Reply to 'Difficulty in the Fermi-Liquid-Based 
Theory for the In-Plane Magnetic Anisotropy in 
Untwinned High-T c Superconductors' 

We show that although our original paper contains 
some unprecise statements there is neither a contradic- 
tion between ARPES and INS in assuming a Sq < 0, nor 
significant physical errors that affect our results. In par- 
ticular, we did never proved Sq > 0; instead, So = —0.03 
had been already used in Fig. 4(b) of Ref. [l|. 

Let us first note, that indeed a positive parameter <5 
would be consistent with the result of a simple quantum- 
chemical calculation yielding t a /tt, ~ (6/a) 4 . On the 
other hand, the LDA Fermi surface (FS) is more com- 
plicated than the one-band model used by us. In par- 
ticular, the saddlepoints are not at (±7r, 0) and (0, ±7r), 
but bifurcated considerably away from this Q. Further- 
more, we have neglected both the bonding plane band 
and the chain band. A recent work by O.K. Andersen 
and co-workers calculating the downfolded plane bands 
reproduce the LDA result in detail, but in a non-trivial 
way: the nearest-neighbor hopping amplitude is smaller 
along the x- than along the y-direction, and the opposite 
is true for longer-ranged hoppings Q. Thus, using the 
notation of Ref. one finds \t\ x \ < \t\ v \. Because the 
first hopping matrix element is most important, we have 
used S < in Fig. 4(b) of Ref. in order to explain 
the anisotropy in the inelastic neutron scattering (INS) 
data by Hinkov et al. [f|. The result |ii x | < \t± y \ makes 
it necessary to re-visit early ARPES data on untwinned 
YBa2Cu306+a; (YBCO) [6]. However, in contrast to our 
Fig. 1 of Ref. 0] in untwinned YBCO, ARPES does not 
observe a closing of the Fermi surface around (0, ±7r). In 
order to see whether a closed or opened Fermi surfaces 
would affect our main conclusions we present in Fig. [T] 
the result obtained from our Fermi-liquid-based approach 
for Sq < and two different tight-binding parameters. 
One clearly sees that, although the FS topology changes, 
the calculated INS response reveals basically the same 
result, i.e. two-dimensional and highly anisotropic with 
two clear maxima along the ^-direction (see also RefQ). 

Finally, if Cooper-pairing is driven by a short-range 
interaction as it is believed in high-T c cuprates, then 
the wave function ip ~ (cosfc x +cosfcj,) corresponds to 
the s-wave component of the superconducting gap func- 
tion. However, as Zhao and Li correctly pointed out [§], 
this admixture of the s-wave component to the original 
d x 2_ y 2-w&ve symmetry does not yield the experimentally 
observed anisotropy of the superconducting gap. In or- 
der to see whether this is a significant physical error we 
present in (e) and (f) the calculated INS response for 
uj = 35mcV. One clearly sees that the results are nearly 
independent of the particular d + s-wave gap representa- 
tions. 

In summary, for So = —0.03 there exists no contra- 
diction between the calculated INS intensity and recent 
ARPES data. Thus, our Fermi-liquid-based theory still 
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FIG. 1: (color online). Comparison of the calculated 
INS intensity (b) and (d) for two different Fermi surface 
topologies (a) and (c) calculated for two different tight- 
binding parameters (tb2, tb3 from Ref. 10]). Different 
types of superconducting gaps are compared in (e) Afc = 
AfJ (cos k x — cos k v ) /2 + Aq (cos k x + cos k y ) /2 and (f) = 
Aq (cos k x - cos k y ) /2 + Ag. 

provides an alternative approach to the stripe scenario in 
order to explain highly anisotropic INS data. 
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